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ABSTRACT 


Computation of Discrete Fourier Transform (DFT) is a challenging task. Especially, on computational 
machines/embedded systems where resources are limited. The importance of Fourier Transform (FT) 
cannot be denied in the field of signal processing. This paper proposes a technique that can compute 
Discrete Fourier Transforms for a matrice or vector with the help of matrix multiplication. Moreover, 
this paper discusses the trivial methods used for computation of DFT along with methods based on 
matrix multiplication used to compute discrete Fourier Transform in addition to the shortcomings. 
The proposed method can help in the calculation of a Discrete Fourier Transformation matrix by 
truncation of values from the proposed generic method which can help in computing DFT of varying 
lengths of vectors. On legacy computing machines and programming environments, having support 
for matrix multiplication, the proposed methodology can be implemented. 
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1 Introduction 


Modern embedded computational systems require flexibility and scalability to cater the abruptly changing demands in 
technology. These demands force the designers to change the system parameters on the fly to ensure the quality of 
service. Fast Fourier transform (FFT) [1] 1s a fast and efficient algorithm for computing the Discrete Fourier transform 
(DFT) [1-{3}. Over the past few years, Image, Audio, Digital communication systems and deep learning-based systems 
have undergone massive production [41124]. Most of the algorithms used for these systems involve the computation of 
Fourier transforms. The computational complexity of DFT 1s a significant factor. Although the Fast Fourier Transform 
algorithm decreased the computational complexity of the Discrete Fourier Transform, the importance of the Discrete 
Fourier transform method still stands in its place. For an educational purpose, a brief report for computing the Discrete 
Fourier transform is presented along with a proposed technique that can compute Discrete Fourier transform using 
matrix multiplications. This paper describes a simple yet effective method for computing DFT for legacy machines and 
devices such as FPGA, portable compute sticks, etc. It is not always feasible to replace legacy computing hardware. In 
some cases where time is not the main concern and functionality is the primary goal, this method can be deployed to 
provide the latest functionality while keeping in view the cos. In other words, this paper provides a software-based 
solution for legacy hardware that enables it to compute DFT with no additional hardware resources. The proposed 
method can be considered as a lookup table. Further, the method for calculating the lookup table, general matrix 
computation, how to compute a DFT matrix, and how to utilize the DFT matrix for computing DFT and Inverse DFT 
is provided as follows. For the computation of DFT, a vector of N elements, a DFT Matrix containing the twiddle 
factors should store V x Nelements. Thus for a vector of MM elements will require to compute a matrix containing 
twiddle factors of size IZ x M 1.e. whenever the size of the input vector is changed a new DFT Transformation Matrix 
needs to be computed every time. Our approach stores a pre-computed general matrix that is similar to a lookup 
table. The number of required elements will be truncated from this general matrix and multiplied with the input vector, 
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which in turn, yields the Fourier Transform for the input Matrix. This technique can be applied in a general computing 
environment and other dedicated machine architectures which support matrix multiplications. 


The main contribution of this paper is that it proposes a method that enables legacy systems to compute DFT using 
the proposed general Matrix scheme. Additionally, some experiments and examples are also presented to validate the 
proposed scheme. A GitHub repository that contains the implementation code along with some real-life examples can 
be found at 


2 Discrete Fourier Transform 


The discrete-time Fourier transform is useful for interconversion of a signal from the spatial domain and frequency 
domain. The frequency-domain has the same information of the signal as the spatial domain but in a different form. A 
discrete Fourier transform of a two-dimensional square matrix can be computed by the following expression. 


The DFT of a vector can be calculated using the following expression: 


N-1 
Xp= do ane PrN EZ (1) 
n=0 
For expressing the inverse DTFT of a vector the following equation can be employed: 
Nai 
So Xpe2™"/N neZ (2) 
k=0 


a 
in = 


The following equation can be used to calculate DFT of a matrix: 


F(u) = BN OS flayew Cfor' u,v € [0,N = 1] (3) 


A two Dimensional matrix can be represented by multiple one dimensional Fourier transforms: 





Pb) = x ear F(z,v)e7"N" for u € [0,N — 1] (4) 


Where Wi = e7J2**/N or e—J27Y/N are the twiddle factors 25}. So, it can be concluded that the DFT of a square 
matrix(// x M) is computed by performing M one dimensional DFTs for each row. Afterward, the same procedure is 
applied to the resultant but column-wise. 


3 Trivial method for calculating DFT 


There are many trivial methods used to compute the Discrete Fourier transform of a vector / Matrix. Equation (1) can 
be used to compute the Discrete Fourier Transform of a given vector x. Typical methods for computing DFT along with 
MATLAB code and simulation are presented. 


The most common method used for computing DFT involves the programming implementation of Equation (1), the 
following program snippet is capable of computing DFT of a provided vector. The code can calculate the DFT of the 
specified NV Values: 


4 DFT Calculation By Equation 


N=15; Define the Number Of points 
x=rand(1,N); % A vector with N Random Values 
DFT=zeros(1,N); 4 Allocation Zeros to a new vector for storing Results 


4 The process (Expansion) 
for K=0:N-1 
for n=0:N-1 
DFT(K+1)= DFT(K+1)+( x(nt+i) * exp((-i*2*pi*K*n) /N)) 
end 
end 


‘https://github.com/abdkhanstd/GSDFT 
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Similarly, the inverse Fourier transform of the vector can also be computed using Equation (2) and the following code: 
AInverse IDFT Calculation By Equation 


N=15; Define the Number Of points 
x=rand(1,N); 4, A vector with N Random Values 
IDFT=zeros(1,N); % Allocation Zeros to a new vector for storing Results 


4 The process Inverse Expansion 
for K=0:N-1 
for n=0:N-1 
IDFT(K+1)= IDFT(K+1)+( x(mt1) * exp((i*2*pi*K*n) /N) ) 
end 
end 
IDFT=IDFT/N; %Now IDFT Contains The inverse 


A DFT transformation matrix can be represented as a Vandermonde matrix as shown below. The elements of this matrix 
are commonly known as twiddle factors: 


a wo oe ao 
ele jel - woe ey 
Pe N N N (5) 
oe go _ oe) 


Where wy = e—N. Whenever a vector is multiplied by the transformation matrix F’. It yields the DFT of the provided 


vector 1.e.: 


A= xX (6) 


The transformation matrix contains NV x N elements, where, N is the total number of elements in the vector k. The 
following MATLAB code can compute the Fourier transformation matrix F’, which, can be used for computing the 
Fourier transformation matrix. 


4Code for Computing DFT Transformation Matrix 


N=8; 4% Transformation Matrix for N Number of Points 
F=zeros(N,N); 
for x=0:N-1 
for y=0:N-1 
F(xt+1,y+1)=(exp((-2*pi*i) /N))~ (x*y) ; 
end 
end 


For computing, the inverse DFT using the transformation Matrix the input matrix X (the DFT of vector x) is multiplied 
to the conjugate of the transformation matrix and by Dividing the resultant with N yields the original signal 1.e. 


1 
— — X fF* 7 
B= (7) 


Another way to compute the transformation matrix is to pre-compute the DFT of an Identity matrix of size IZ x M where 
M is the size of the vector x whose DFT needs to be calculated, the resulting matrix can be used as a transformation 
matrix itself. The following expression explained the calculation of the transformation matrix. 


F = DFT (Imxm) (3) 


This method will yield the same transformation matrix as Vandermonde matrix explained above section. For computation 
of DFT for a two dimensional matrix typically the Fourier transform can be computed by computing the fourier transform 
of each column of the matrix x, then again computing the Fourier transform of each column of the transpose of the 
resultant and then again transposing. Mathematically it can be written as: 
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X = DFT (DFT(z)') (9) 


If it is desired to calculate the DFT using the transformation matrix, no transposing 1s required the DFT can be calculated 
by simple matrix multiplication of the input matrix with the transformation matrix which can be represented as: 


AS 6g XT (10) 


For computing, the inverse DFT using the transformation Matrix the input matrix X (the DFT of vector x) is multiplied 
to the conjugate of the transformation matrix and by Dividing the resultant with N yields the original signal 1.e. 


1 
a2, i 11 
NT (11) 
Another way to compute the transformation matrix is to pre-compute the DFT of an Identity matrix of size M x M where 
M is the size of the vector x whose DFT needs to be calculated, the resulting matrix can be used as a transformation 


matrix itself. The following expression explained the calculation of the transformation matrix. 


F = DFT (Imm) oa 


This method will yield the same transformation matrix as Vandermonde matrix explained above section. For computation 
of DFT for a two-dimensional matrix typically the Fourier transform can be computed by computing the Fourier 
transform of each column of the matrix x, then again computing the Fourier transform of each column of the transpose 
of the resultant and then again transposing. Mathematically it can be written as: Where X is the resulting two- 
dimensional DFT, F’ is the transformation matrix and z is the input square matrix. 


4 Proposed Method 


A problem arises with the transformation matrix that DFT can only be computed for the size of NV points. If it is required 
to compute DFT of points other than VV, The transformation Matrix has to be computed again, the DFT computation 
is already compute-intensive. Computing another transformation matrix for this new-sized vector will require more 
computation and storage in the memory. 


There is a possibility, that a general DFT Transformation matrix can be employed like a lookup table (which is larger). 
To construct a generic matrix, a mechanism is required for truncating values from a bigger pre-computed matrix. 
Suppose that the Generic Transformation matrix is computed for a VX N. This matrix will be capable of computing 
the DFT of N point vector. Now a situation occurs that it is required to compute the DFT of V/ Points vector Where 
M < N asimple truncation can be done in the generic matrix already computed for NV points. A new matrix can be 
produced from a larger generic matrix, the truncation can be done from the first row till 1/*” row and from the first 
column till M‘” column. 





NxN 


Figure 1: Truncation of pre-calculated values from larger Generic Matrix for computation of DFT. 


In Figure[I] the bigger bold Square represents the generic matrix, and the dotted square represents the smaller truncated 
new transformation matrix. The purpose of using this mechanism is to provide a flexible and scalable stored table DTFT 
calculation Matrix. A twiddle factor is stored in the transformation matrix. Here, we can see that it holds the effect 
of N. The goal is to create a general DTFT matrix, that can be further truncated according to the required points. To 
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achieve this, we have to remove the effect of V in the equation. Instead of storing the value, we can simply store the 
value —27r2 after truncating. The new matrix can be divided by the total number of new points M/ and afterward taking 
the exponential. This yields a new transformation matrix. 


Trivially, the twiddle factor is represented as: 





Where: . 
a ee (14) 
By taking Log on both sides: 
log(w) = log Ca (15) 
Simplifying Equation[1I5]yields: 
—271 
l = 16 
og(w) = (16) 
Nlog(w) = —27i = Wnew = —6.2832% (17) 


Here, instead of using the twiddle factor, we simply store —272 instead of the twiddle factor in the general DFT 
transformational Matrix. For our transformation matrix, we store the new twiddle factor wy-, . Then the generic 
Transformation Matrix can be represented as: 


Wnew -0.0 Gane a0 bee Whew -0.(N — 1) 
Deu lO) one oe — Whew -1.(N — 1) 
G= ; | (18) 
whew -(N—1).0 whew -(N—1).1 22. thew *(N—1)-(N-1) 


G represents the new generic transformation matrix, this matrix does not have the effect of the total number of elements, 
thus can be truncated for different sized matrices. After truncation, some more processing is required. The new truncated 
transformation matrix needs to be divided by NV, which, is the total number of elements for the vector whose DFT is to 
be calculated: 


F=ewn (19) 


The point to be noted is that the number of operations is almost the same as required by the computation of the regular 
transformation matrix. The above procedure can mathematically be explained by the following set of equations. 


Nlog(w)  —2zi 











= 20 
N N (20) 
_ On 4 
N log(w) = ——~ (21) 
Taking the Exponential of the above equation the expression reverts. 
G=2 (22) 


This generic table holds the values that are the exponents. A Matlab simulation for the proposed method is provided to 
support our claim: 


ASimulation: DFT of vector using Matrix Multiplication 
Generating a general DFT Matrix that can be cropped according to a new 
General Matrix row and columns 


N=10; 

v=rand(1,N); % V is the vector 

Gr=256; 4, Dimension of the Generic Matrix 
Gc=Gr; 


% Calculation The Generic Matrix 
for x=0:Gr-1 
for y=0:Gc-1 
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G(x+1, y+1)=(-6.2832+*i*x*y) ; 
end 
end 


/ Truncating Transformation matrix from Generic matrix 
F=exp(G(1:N,1:N)/N) ; 

4» Calewlating DFT 

DFT=v+F 


4.1 Experiments and results 


Likewise a traditional research paper, the proposed method does not require traditional experiments. However, to 
validate the proposed method, some results are computed on random vector and 2D vector (matrices). The results show 
that the proposed method is fully capable of computing DFT. 


Some simulation and experimental results for small vectors are presented in Figure|2| In Figure 2, the first column 
shows the original signal. The second column represents the frequencies and amplitude obtained using the standard 
DTFT method. Moreover, the third column shows the results obtained by the proposed DFT calculation method. It can 
be seen that the obtained results are the same as the trivial DTFT calculation method. A code for this example can be 
found in the GitHub repository as Example 1, Example 2, and Example 3 for rows 1, 2, and 3 respectively. The first row 
presents the results for a signal with two fundamental frequencies 1.e., 50 and 120 Hz with 150-time samples. While, 
the second row presents the results for a signal with five fundamental frequencies 1.e., 10, 50, 120, 150, and 180 Hz 
with 150-time samples. Moreover, the third row presents the results for a signal with five fundamental frequencies 1.e., 
10, 50, 120, 150, and 180 Hz with 1500-time samples. 
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Figure 2: The simulation results of a sinusoidal signal with basic frequencies between 10 and 180Hz. 
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Figure 3: Some selected samples of DTT computation of images using the proposed method and the trivial DFT 
computation method. 


The results for the computation of 2D DFT are presented in Figure [3] In Figure 3, the first column represents the 
original images, the second column represents the 2D DFT computed by trivial 2D DFTcomputation method and, the 
third column represents the 2D DFT computed using the proposed method. The first image in Figure 3-(a) is a black 
and white image with a size of 30 x 30 pixels only. Moreover, in Figure 3-(d) a repeated pattern sample is shown with 
an actual size of 200 x 200 pixels. Further in, Figure 3-(e) and Figure 3-(f) the shifted frequency maps are shown which 
were obtained using trivial and the proposed DFT computation method for 2D vectors. An image of size 50 x 50 is 
presented in Figure 3-(g). Whereas, Figure 3-(h) and Figure 3-(1) the shifted frequency maps are shown which were 
obtained using trivial and the proposed DFT computation method for 2D vectors. From the results, it can be seen that 
the proposed method is capable of producing similar results while using lesser hardware resources. 


4.2 
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Conclusion 


The method proposed in this paper can help in the calculation of a Discrete Fourier Transformation matrix by truncation 
of values from the proposed generic method which can help in computing DFT of varying lengths of vectors. The 
proposed method can be considered as a lookup table that enables the computiation of DFT. The method and the relevant 
code are provided in the manuscript. Moreover, the results show that the proposed method is fully capable of computing 
DFT. On legacy computing machines and programming environments, having support for matrix multiplication, the 
proposed methodology can be implemented. 
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